Proposal for generating and detecting multi-qubit GHZ states in circuit QED 
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We propose methods for the preparation and entanglement detection of multi-qubit GHZ states 
in circuit quantum electrodynamics. Using quantum trajectory simulations appropriate for the 
situation of a weak continuous measurement, we show that the joint dispersive readout of several 
qubits can be utilized for the probabilistic production of high-fidelity GHZ states. When employing 
a nonlinear filter on the recorded homodyne signal, the selected states are found to exhibit values 
of the Bell-Mermin operator exceeding 2 under realistic conditions. We discuss the potential of the 
dispersive readout to demonstrate a violation of the Mermin bound, and present a measurement 
scheme avoiding the necessity for full detector tomography. 
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I. INTRODUCTION 

Performing measurements on separable input states is 
an elegant way to create entangled states of two or more 
qubits without the need for high-fidelity two-qubit gates 
[J, 0, 0, 0, Si- For this method to work, the measure- 
ment needs to be a joint readout of several qubits, as 
single-qubit measurements are insufficient for entangle- 
ment generation. In circuit quantum electrodynamics 
(cQED) 01, the dispersive readout constitutes such 
a multi-qubit measurement. The qubits, realized by su- 
perconducting charge qubits, are coherently coupled to 
the voltage inside the resonator and give rise to a state- 
dependent shift of the resonator frequency. Heterodyne 
or homodyne detection of the phase of microwave radia- 
tion transmitted through the resonator can thus be used 
to infer the state of the qubits. Ideally, the correspond- 
ing measurement operator is simply a weighted sum of 
the qubits' Pauli operators cr^ , where the weights are 
conveniently adjusted by the detunings of the respective 
qubits from the resonator frequency. 

The idea of probabilistic state-preparation by measure- 
ment has recently been applied to a 2-qubit cQED system 
by Hutchison et al. [1]. Their theoretical study included 
the adverse effects of qubit relaxation and dephasing, and 
has shown the practical applicability of the method, even 
for realistic decay and decoherence rates as currently 
realized in cQED experiments. Here, we extend this 
method to the generation of multi-qubit Greenberger- 
Horne-Zeilinger (GHZ) states @. For superconducting 
qubit systems there have been successful demonstrations 
of Bell-state preparation [l^, [ill, j ^^'^ proposals for 
creating GHZ states, mainly focusing on the generation 
via two-qubit gates and qubit-qubit interactions, see e.g. 
[H, . Instead of employing such entangling gates for 
generating a GHZ state, we propose a scheme tailored to 



cQED, consisting of one-qubit rotations and a dispersive 
measurement only. Based on quantum trajectory simula- 
tions, we show that currently attainable values for qubit 
decoherence and decay allow for the creation of three- 
qubit GHZ states in cQED with high fidelity and high 
degree of entanglement. The degree of entanglement can 
be increased at the cost of lowering production rates. 

In order to verify the production of the desired GHZ 
state, we propose to use a second dispersive readout. Be- 
cause the GHZ state is maximally entangled, this verifi- 
cation is related to proving the violation of a Bell-type 
inequality [l^ [13 ■ However, proving such a violation in 
a loophole-free fashion turns out to be a much more chal- 
lenging task in cQED. Given the required measurement 
time of hundreds of nanoseconds, space- like distances (in 
the sense of special relativity) are of the order of tens 
of meters and thus difficult to achieve in a cQED setup, 
and the dispersive readout is in fact inherently nonlo- 
cal. Accepting that the locality loophole therefore can- 
not strictly be closed, we discuss the potential of the 
dispersive readout for observing quantum correlations in 
a 3-qubit GHZ state, as well as the potential for devis- 
ing a factorizing measurement that is local in the no- 
signalling sense [40| . Using quantum trajectory simu- 
lations including the measurement imperfections caused 
by qubit decay, we show that a convincing violation of 
the Bell inequality would require a signal-to-noise ratio 
which is currently out of experimental reach, but may be 
approached once efficient methods for protecting qubits 
from decay have been devised, or with improvements in 
the noise performance of microwave amplifiers 18]. 

Our paper is organized as follows: In the following 
section, we present the central idea of generating and 
detecting multi-qubit GHZ states by dispersive measure- 
ments, starting with the exposition of our proposal in 
the idealized situation of no qubit decoherence and de- 
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cay. The rest of the paper is devoted to the consequences 
of measurement imperfections introduced by decay dur- 
ing the measurement process. In Section IIIH we specify 
the treatment of qubit decay and continuous homodyne 
detection using an effective stochastic master equation 
previously introduced in Ref. [1]. Quantitative results 
from solving this master equation for the situation of 
GHZ-state generation are presented in Section IIVI We 
describe different protocols for accepting or rejecting a 
generated state as a GHZ state, and show in particular 
that nonlinear filtering offers a significant advantage over 
simple boxcar filters. In Section |V] we discuss the detec- 
tion of GHZ states within the dispersive measurement 
scheme and comment on the potential to violate a Bell- 
type inequality. Finally, our conclusions are presented in 
Section ED 



modyne signal for the microwaves transmitted through 
the resonator. In the absence of qubit decay and de- 
coherence, the probability distribution p{s) for the in- 
tegrated signal s (see Sec. IIIII for precise definition) 
takes the form of Gaussian peaks, which initially over- 
lap stron gly and separate with increasing measurement 
time t m. Ha, ill. In the limit of negligible over- 

lap between peaks, the dispersive readout corresponds to 
a projective measurement of the operator A — J2j ^j'^j^ 
where the weights 5j = Xj/x the fractional contri- 
butions to the mean dispersive shift, x — '^j Xj ■ In 
detail, the preparation scheme can be described as fol- 
lows: 

(i) Arrange all qubit dctunings such that the system 
is dispersive, and mutual qubit dctunings are large 



II. IDEALIZED PREPARATION AND 
DETECTION OF GHZ STATES 

In this section, we lay out the essential ideas behind the 
preparation and detection of GHZ states using the joint 
dispersive readout typical for cQED. To keep the discus- 
sion as clear as possible, our exposition in this section will 
ignore the adverse effect of qubit decay and decoherence, 
and other possible sources of measurement imperfections. 
We turn to the full discussion of the realistic situation in- 
cluding these effects in the following sections. 

The GHZ state Q is the maximally entangled multi- 
qubit state of the form 
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where | ■)j denotes the state of qubit number j. GHZ 
states have received much attention in the context of vi- 
olation of Bell-type inequalities, see e.g. [l^, [13, HH IHj 
[23} . ruling out classical local- hidden- variable theories as 
a valid description of nature. 

Partly, the beauty of the GHZ state lies in the fact 
that, in principle, violation of classicality can be proven 
with a single measurement of the corresponding Bell- 
Mermin operator M [l^ , see Eq. (U) . This has to be con- 
trasted with the situation of the two-qubit CHSH scheme 
[U [25! , where such a proof necessarily requires accumu- 
lating statistics. Key to this difference is the property of 
the GHZ state being an eigenstate not only of M, but 
also simultaneously of the measurable parity operators 
which sum up to M. The 7V-qubit GHZ state is known 
to violate a Bell-type inequality by an amount that grows 
exponentially in the number of qubits [l9| . 



A. Preparation scheme 

In the dispersive readout of cQED, the measurement 
outcomes are inferred from the detection of the ho- 
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FIG. 1: (Color online) Dispersive measurements employed 
for (a) generating a GHZ state, and (b) detecting the par- 
ity n = Ylj <^j- Both panels show the probability density 
p{s) for the integrated homodyne signal s for the concrete 
example of a 3-qubit system, (a) For the generation of a 3- 
qubit pre-GHZ state, the dispersive shifts are fixed at ratios 
Xi • X2 : X3 = 1 ■ 1 • 2. Ideally, the Gaussian peaks belong- 
ing to the 5 measurement results {±4, ±2, 0} separate with 
increasing measurement time t (here: t = S/Fci), allowing 
for a rehable projective measurement when using appropri- 
ate thresholds, e.g. ui, 1/2 for the selection of the measure- 
ment outcome "0". (b) The dispersive parity measurement 
requires identical dispersive shifts, Xi • X2 : X3 — 1 • 1 • 1- 
The four measurement outcomes a,; G {±3, ±1} then allow 
the inference of the parity value by Hi — — sin(ai7r/2). 
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compared to the qubit-qubit interaction strengths. 
The initial state is the state with each qubit in its 
ground state, 

(ii) Perform 7r/2 rotations on each of the N qubits, 
preparing the state 2-^/^ (g)^^^(| i + 1 1 ),)• 

(iii) Keeping the system dispersive and mutual qubit 
detunings sufficiently large, adjust the qubit de- 
tunings such that their dispersive shifts assume the 
ratio, 

Xi ■X2- ■■■■■ XN~i : XN ^ I ■■ I ■■ ■ ■ ■ ■■ I ■■ N ~1, (2) 

and perform a dispersive measurement. Ideally, 
this corresponds to a projective measurement of the 
observable A = J^j ^j'^j- Conditioned on the mea- 
surement result being "0", see Fig. [IJa), we thus 
obtain the pre-GHZ state 

IpGHZ) = (|ii...iT> + ITT...Ti))/A/2. (3) 

(iv) In the final step, a tt rotation is applied to qubit 
A^, yielding the GHZ state, Eq. Alternatively, 
one may choose a different computational basis by 
interchanging the "t" and "J," labels for qubit N. 

The necessary adjustment of the Xj ratios is pos sible in 
cQED samples employing local flux-bias lines flSl, which 
allow for the fine tuning of individual qubit frequencies. 
We note that the scheme requires the resolution of only 
~ 2N different peaks, which should be compared to the 
need for application of {N — 1) two-qubit gates for the 
preparation of the same GHZ state via gates, see e.g. 

B. Detection scheme 

Ideally, the confirmation of the GHZ state production 
and the verification of its quantum correlations proceed 
by a measurement of the Bell-Mermin operator 

, N N s 

M^2--H(Yla--YlaA. (4) 

For the A^-qubit GHZ state, this operator takes on the 
value 2^~^, while local-hidden variable theories predict 
an outcome < 2^/^ if N is even, and < 2(^-i)/2 if TV is 
odd [l^ - thus leading to a violation that grows expo- 
nentially in the qubit number. 

In the general case, the Bell-Mermin operator is not 
amenable to a direct measurement. However, for N 
qubits, it can be decomposed into 2^~^ 7V-qubit par- 
ity operators, which are more easily accessible by experi- 
ment, and the GHZ state is a simultaneous eigenstate of 
all the relevant parity operators. The specific form of the 
Bell-Mermin operator in the three-qubit case is given by 

M = gIgIgI - alalol - afa^a^ - afa^a^ (5) 



In the ideal case, one would perform the 2^^^ parity 
measurements, using a quantum non-demolition method 
on one and the same state and not requiring repeated 
measurements. 

Since the dispersive readout does not realize exact par- 
ity measurements, we will accept the necessity to repeat 
measurements and acquire statistics. Instead of the par- 
ity, the dispersive readout can easily access the operator 
A = J2j '^j I from which the value of the parity Ylj can 
be uniquely inferred, see Fig.[ljb). Using single-qubit ro- 
tations mapping the appropriate x and y axes to z (2^ . 
all the required parities can be measured dispersively. 

The crucial step thus consists in tuning all dispersive 
shifts to be identical. As before, this can be achieved 
by adjusting qubit detunings using local flux-bias lines. 
Compared to the setting employed for the GHZ state 
generation, it is in fact only the detuning of the A^- 
th qubit that needs to be changed. Ideally, the mea- 
surement of A then leads to the measurement outcomes 
a, G {±N, ±{N -2),..., ±e}, terminating with ^ = 1 if 
N is odd and with £ = if is even. The inferred 
parity outcomes simply alternate in sign according to 
Hi — ~ sin(ai7r/2) for odd A^ and by H^ — cos(ai7r/2) for 
even A^. It is important to note that, while the number 
of required different measurements grows exponentially 
with A^, the number of measurement outcomes that need 
to be resolved is given by A^ -I- 1, only growing linearly 
with the qubit number. This should be compared to 
the situation of a full state readout, which would require 
resolution of 2^ different peaks and dispersive shifts to 
be spread over an exponentially large frequency range, 
Xj =2^X0- 

Both the generation and detection scheme will obvi- 
ously suffer from qubit decoherence and decay. The sub- 
sequent sections take into account these effects and study 
quantitatively how the idealized proposal performs under 
more realistic conditions. 



III. MODEL 

For the generation and subsequent detection of a multi- 
qubit GHZ state we consider a cQED system 0] con- 
sisting of three superconducting charge qubits coupled 
to the fundamental mode of a microwave resonator. The 
model of the system and notation follow those in Ref- 
erence Neglecting the possible influence of levels be- 
yond the two-level approximation for the superconduct- 
ing qubits, the system is described by a driven Tavis- 
Cummings Hamiltonian (30j 

3 J 

(ae*e''^"* -fa^ee-'""*), (6) 

where we set h = 1, uj^l2'n denotes the resonator fre- 
quency, and e the strength of the measurement drive. 
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The qubit frequencies u}q,j/2TT are considered to be tun- 
able individually, as realized by local flux-bias lines in 
recent cQED experiments . The qubit-resonator cou- 
plings are given by gj , whose signs are determined by the 
location of the respective qubit within the resonator. For 
concreteness, we will focus on the case of a A/2 coplanar 
waveguide resonator, with two qubits placed close to one 
end, and the third qubit on the opposite end, leading to 
a relative sign sgn(5i) = sgn(g2) = - sgn(53). 

The system is to be operated in the dispersive regime, 
where = — cj^l = \9j/^i\ 1- Under these 

conditions the interaction term in Eq. ^ can be adiabat- 
ically eliminated @ , such that the effective Hamiltonian 
in the frame rotating with the measurement drive fre- 
quency bJm reads 



3 



3 



(7) 

where the detuning between measure- 

ment drive and resonator, and y,- = 9^ / denotes the 
dispersive shift due to qubit j [4l| . Here, the qubit-qubit 
coupling ~ J via virtual photons has been neglected, 
as is appropriate for sufhcient detuning between qubits, 
J ^ — Aj/|. The effects of qubit decay and cavity 
photon leakage are taken into account within a master 
equation description. Specifically, we include intrinsic 
qubit relaxation with rates 7ij, Purcell-induced relax- 
ation with rates ^pj [sH . Is^ l , and photon decay from the 
cavity with rate k. Pure dephasing can be strongly sup- 
pressed by proper design of the superconducting qubit 
[331, and will be neglected here. (We have checked that 
inclusion of pure dephasing at small rates, comparable 
to those achieved in [33|, does not significantly alter our 
results.) 

As demonstrated in Ref. 8, one can dramatically sim- 
plify the resonator-qubit master equation and reach an 
effective master equation for the qubits only, given that 
photon decay is fast. Specifically, we assume that A^ = 0, 
and require 



K> max{e,^^. |xj |} 



(8) 



Under these conditions, an analogous separation of qubit 
and resonator degrees of freedom can also be reached on 
the level of the stochastic master equation (SME), ap- 
propriate for the situation of continuous homodyne de- 
tection of the emitted microwave radiation Q. The ef- 
fective SME for the qubit density matrix pj conditioned 
on the measurement record 



(9) 



is given by 



PJ = Cpj + v^C(t)Al 5,cj'^]pj, (10) 

3 

where we are using notation identical to Ref. H: M[c\ 
is the measurement operator given by M[c\pj = (c — 



(c))pj/2 + pj{c— (c))/2, ^(t) represents Gaussian white 
noise with zero mean and {£,{t)£,{t')) = S{t—t'), and Yd — 
r/Fin denotes the effective measurement rate, reduced by 
an efficiency factor _with respect to the maximum rate 
64x2 1 



r 



42l | . The Liouvillian is defined as 



Cp 



^3. 



4xkl 



J2s,<j'-,p 



(11) 



with the measurement-induced dephasing rate Fd = 
Fm/2, and the usual dissipation superoperator 'D[A]p = 
ApA'' — {A^ A, p} /2. Assuming mutually distinct qubit 
frequencies, we have treated the Purcell effect in the sec- 
ular approximation. Specifically, we neglect the cross- 
terms in 2? [erf -l- — cr^] which are interference effects 
for radiation from different qubits, and which only be- 
come important if the qubit frequencies are sufSciently 
close, i.e. \Ai — Aj\ <^ ^pi^ 



|34| . In our case, Purcell- 
induced decay and intrinsic decay can be treated on the 
same footing, and in the following we will assume similar 
decay rates for all qubits and subsume them under the 
shorthand 7 — jij + jpj . Finally, the integrated signal s 
is simply given as the time integral of the measurement 
record for the full measurement time t, s = Jg dt' J{t'). 



IV. PREPARATION OF THE GHZ STATE 
UNDER REALISTIC CONDITIONS 



We now turn to the situation of GHZ state prepara- 
tion in the presence of qubit decay, which we study using 
quantum trajectory simulations based on the stochastic 
master equation (|10p . Following the steps (i)-(iii) de- 
scribed in Section |TT1 the system is initialized and dis- 
persive shifts are adjusted for the measurement step. 
The interplay of measurement-induced dephasing, grad- 
ual state projection, and the simultaneous qubit decay 
are captured by the conditional density matrix pj, where 
each simulation run generates a particular measurement 
record J{t) up to a final measurement time t, correspond- 
ing to the experimentally accessible homodyne signal. 

Since preparation of the correct pre-GHZ state is prob- 
abilistic (ideally, state generation succeeds with proba- 
bility P ~ 1/4 in the present case), one has to define a 
criterion ("filter") for success of preparation, and post- 
select the corresponding subensemble |2^. In principle, 
the information available to the filter is the full measure- 
ment record. In the following, we will discuss two differ- 
ent filters, the linear boxcar filter and the full nonlinear 
Bayesian filter and compare their performance in select- 
ing high-fidelity GHZ states under realistic conditions. 

The simple filter already outlined in Section |TT] is 
the linear boxcar filter. It compresses each measure- 
ment record into a single number, the integrated signal 
s = dt' J{t'), and declares successful pre-GHZ state 
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FIG. 2: (Color online) (a) Histogram of the integrated signal 
after a measurement time t = S/Fci, and probability distri- 
bution p{s) in the absence of any decay (red/gray curve), (b) 
Scatterplot (blue/gray dots) showing the correlation between 
the expectation value of the Mermin operator (M) and the 
integrated signal s for t = S/Fci. Each point corresponds to 
one of 10 000 trajectories. For comparison, the correlation in 
the ideal case of no decay is shown as the red/gray curve. 
The boxes indicate the action of the boxcar and the nonlin- 
ear filtering scheme. Parameters are chosen as Fd = Fci/2, 
7/Fci = 1/35 for j = 1, 2, 3, and (Si = (S2 = 3/4 and 5^ = 3/2. 



preparation whenever s falls within the limits of appro- 
priately chosen thresholds, 1^1 < s < 1/2- Otherwise, the 
state is rejected. 

The results for the integrated signal of many such 
measurements are conveniently plotted in form of a his- 
togram, see Fig. [21a). When compared to the probability 
distribution expected in the ideal case of no decay, one 
observes that qubit decay leads to a distortion of the 
probability density with an overall shift of probability 
density towards the left-most peak, i.e., towards the sig- 
nal associated with the ground state. The shift is thus 
easily understood as a consequence of decay processes 
acting during the finite measurement time. 

As a benchmark for the quality of the generated states 
and its correlation with the integrated signal, Fig. [D^b) 
shows a scatterplot of the expectation value of the Bell- 
Mermin operator (M) versus the integrated signal for 
10 000 individual measurement trajectories. For compar- 
ison, the corresponding scatterplot in the ideal case of no 
decay is shown to collapse to a single curve. The scatter 
in the nonideal case results in trajectories of the same 
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FIG. 3: (Color online) Time traces of the signal J(t) for in- 
dividual quantum trajectories. The traces are smoothed over 
time O.IF^^ (cyan/light gray) and F^^ (blue/dark gray). For 
(a) and (b) the expectation of the Mermin operator is large, 
{M) > 3.9, whereas for (c) and (d) it is small, (M) < 0.1. The 
horizontal lines indicate the values J{t) would take on average 
for the integrated signal s to be at the peaks of Fig.[2ja). All 
4 traces are selected by boxcar filter on the integrated signal, 
such that they all lie close to the middle of the center peak. 
For (b), (d) the relaxation is low, Fci/7 — 142, and trajecto- 
ries with extremal values of (M) can be distinguished by eye. 
For (a), (c) the measurement time is shorter and relaxation is 
faster Fci/7 = 35, nevertheless the nonlinear filter is still able 
to reliably estimate (A/) , as is demonstrated in Fig. 3] 



integrated signal, but very different values of (M), and 
thus in a significant number of falsely accepted states 
within the simple boxcar filtering. 

The essential mechanism for false acceptance of states 
is illustrated in Fig. [3l showing the measurement record 
J(t) as a function of time for four individual trajectories. 
Very roughly, the trajectories with integrated signal s 
close to can be divided into two categories: trajecto- 
ries with measurement records J{t) fluctuating around 
J{t) — 0, see Fig. [3Ia),(b) and measurement records 
showing larger variations of J{t) which accidentally av- 
erage to s = upon integration. Trajectories of the first 
category correspond to the correct pre-GHZ state with 
high probability. On the other hand, an example from 
the second category consists of trajectories which, with 
high probability, initially assume the state | itt ) with 
(A) = 2, and then suffer a decay process in qubit 3 at 
some intermediate time, thus transitioning to the state 
UTi) with (A) = -2, see Fig. [lc),(d). 

This insight also points to a remedy for the boxcar 
filter. The full measurement record can, when spaced 
densely enough, be used to reconstruct the actual un- 
derlying quantum trajectory pj(t) in the following way: 
Given that the state before the onset of the measurement 
[see step (ii) in Sec.|ll] as well as the parameters entering 
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FIG. 4: (Color online) Expectation value of the Mermin op- 
erator (M) as a function of acceptance probability, for sev- 
eral ratios 7/rci; parameters are chosen as in Fig. 2. Solid 
(dashed) lines show the results using the nonlinear (boxcar) 
filter. (See text for details.) Using nonlinear filtering, the 
fraction of accepted trajectories with high {M)-vaIue can be 
substantially increased. For an acceptance probability < 1/4 
the advantage of the nonlinear scheme becomes apparent. For 
each point, (A/) is obtained by averaging over 20 000 trajecto- 
ries and optimizing with respect to measurement time t and 
boxcar thresholds. The inset shows the expectation value (M) 
as a function of the ratio rci/7 for acceptance probability 
of 1%. 



the stochastic master equation are known with sufficient 
accuracy, one can successively determine the Wiener in- 
crements dW{t) = ^(t)dt from the measurement record. 
These, in turn, can then be used to propagate pj from the 
initial time to the measurement time t, and the result- 
ing Pj (t) encodes the expected value of the Bell- Mermin 
operator via (M) — tr[pj{t)M]. This procedure corre- 
sponds to a nonlinear filter [2a |. with an acceptance cri- 
terion based on the value of (M) itself, see Fig. [^fb). 

The advantage of using the nonlinear filter is high- 
lighted by Fig. 31 which compares the performances of 
boxcar and nonlinear filter. For acceptance probabil- 
ities smaller than the ideally attainable P = 1/4, we 
find that the nonlinear filter constitutes a significant im- 
provement over the boxcar filter. Specifically, for ratios 
rci/7 ^ 4 currently supported by experiments, the non- 
linear filter will be crucial in order to reliably exceed the 
value (Af) — 2, which is the relevant Mermin bound for 
violation of local-hidden variable theories in this case. 
Figure [4] demonstrates that, when exploiting the trade- 
off between large expectation values of {M) and high ac- 
ceptance probabilities, high-fidelity GHZ states can be 
prepared under realistic conditions. 



V. GHZ STATE DETECTION UNDER 
REALISTIC CONDITIONS 

The measurement of the Bell-Mermin operator via par- 
ity detection, presented in Section |TT1 requires the resolu- 
tion of ~ peaks in the probability density p{s) of the 
integrated signal. While clearly advantageous relative to 
the resolution of ~ 2^ peaks needed for a full readout, 
the parity detection remains difficult with current experi- 
mental parameters due to the qubit relaxation within the 
measurement time. In the following, we discuss a scheme 
that avoids this problem. 

The key of this scheme lies in the fact that at low 
temperatures, decay into the state | ff) = | TT ' ' ' T ) is 
negligible. We note that this is similar to Kofman and 
Korotkov's use of the "negative result outcomes" to avoid 
the effects of measurement crosstalk in Bell tests using 
superconducting phase qubits ^2§\. False positive events 
in the detection of the state | ) can thus be suppressed 
by setting the acceptance threshold v for the integrated 
homodyne signal sufficiently high. Using this insight, we 
construct a measurement B by assigning the measure- 
ment outcomes "0" , "1" to the cases where the signal is 
respectively smaller or larger than a preset threshold. We 
can describe B in the language of generalized observables 
as a positive operator valued mapping (POVM) jsH], by 
specifying its effects 

Ei=a\t){il\, (12) 
En = l-Ei. (13) 

Here, a = P\^j{s > v) is the probability that the signal 
exceeds the threshold v given the system was prepared in 
I ) . This probability is set by the decay of the | fl" ) state 
during the measurement time, and is analogous to the 
detector efficiency in quantum optics. As a result, 1 — a 
can be described as a "false negative" probability that 
the measurement fails to detect a valid | i[ ) state. Exper- 
imentally, a can be determined by repeatedly preparing 
the system in Iff) (using single-qubit tt rotations), and 
subsequently performing the measurement. This proce- 
dure yields the expectation value \B\ -ff), which is 
identical to the fraction of the cases where s > ly, and 
hence to a. In general, complete characterization of a 
POVM via detector tomography [sl, [13, [M., JM] requires 
many measurements and a numerical optimization pro- 
cedure to ensure the resulting POVM remains physical. 
Due to the simple structure of the measurement B, it 
may be conveniently characterized by determining only 
a single parameter a. 

The measurement B can now be combined with single- 
qubit rotations to determine the parity. We perform all 
combinations of n-qubit bit flips, < n < A'^, and sum 
the measured (B) with relative sign (—1)". For clarity 
we specialize to the 3-qubit case, and define 

-(B) - Ki?f7?) - {a^.Ba^) - {alBal) 

+ (aW^Ba^ai) + {ala^Bala^) + (aj^af Ba^af ) 
- (arafafSarafag-). (14) 
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FIG. 5: (Color online) False negative probability 1 — a and 
worst-case value for the false positive probability /3 (definition 
see text) versus threshold u. The horizontal line indicates the 
necessary constraint on a to violate the Mermin inequality. 
The measurement time is chosen as t — S/Fci. 

The value of fzzz is proportional to the parity measured 
in the z-basis, fzzz — Q;(o'iO'|crf), with the proportional- 
ity constant being a as defined above. It is straightfor- 
ward to extend this scheme to the actual parities required 
for determining the value of the Bell-Mermin operator by 
prepending additional single-qubit rotations. 

The expectation of the Bell-Mermin operator can now 
be related to the actual measurements via F — a{M), 
where 

F — fxxx fxyy fyxy fyyx- (-i^) 

Thus, the measurement of the 32 expectation values en- 
tering into F and determination of a allow for the extrac- 
tion of (M) = F/a, with no restrictions on the qubits' 
decay rates [43| . 

As explained in Section |T] the nature of the dispersive 
measurement prevents us in principle from a strict vi- 
olation of a Bell-type inequality. However, in the limit 
where the measurement effects factorize into tensor prod- 
ucts over the single-qubit Hilbert spaces, i.e. 

E,,k = E^^^ ® Ef ® E^^\ (16) 

the measurement can be considered local in the sense 
of the no-signalling property [13] ■ The effect defined in 
Eq. (fT^ obeys such a factorization 

^^i-"(IT)i(Tli)® (IT)2(Tl2)®(IT>3(tl3)- (17) 

Similarly, the rotated measurements entering into F fac- 
torize in this sense, provided the rotations themselves 
also factorize. This additional requirement holds not only 
for perfect single-qubit rotations 29] , but also for imper- 
fect rotations, as long as there is no coupling or crosstalk 
between qubits during the rotation pulse. For example, 
independent single-qubit relaxation processes during a 
finite-duration rotation pulse do not spoil the factoriza- 
tion property. By contrast, a rotation of qubit b caused 



by a rotation pulse on qubit a no longer factorizes. In the 
following, we will assume that such crosstalk is negligi- 
ble. In that case, the argument of Mermin applies, which 
states that a local hidden variable theory has bounds on 
the allowed F, —2 < F < 2 ^1^. Meanwhile quantum 
mechanics allows for (Af) — 4 and hence if a > 1/2 there 
is the possibility to violate Mermin's version of the Bell 
inequality. 

Figure O shows the variation of the false negative 
probability (1 — a) with threshold v, so that the re- 
quired threshold for a > 1/2 can be read off. Since 
the derivation of the Bell inequality required factoriza- 
tion of measurement effects, we estimate the corrections 
to Eq. p^ . In our case, the largest correction will 
be due to misidentification of states from the subspace 
{| ITT ) , I rn ) , I Tit >}, for which A = a| = 1. We 
put an upper bound on this misidentification probability 
(3 by assuming that there is no decay out of this sub- 
space and thus assume that Pa=i{s), the distribution of 
the homodyne signal arising from this subspace, is Gaus- 
sian. Under these conditions, one obtains a worst-case 
estimate of the "false positive" probability /3 as a func- 
tion of v. 

Figure [H] shows that with a low rate of qubit decay, 
rci/7 = 20, we find a > 1/2 and a low probability of false 
positives, (3 ~ 0.002, meaning that a meaningful violation 
of a Bell-type inequality should be possible. Conversely, 
for a more realistic rate of qubit decay rci/7 = 5, the 
requirement a > 1/2 leads to significant false positive 
rates (3 ~ 0.16, and factorization of Ei breaks down. 
We note that the required rci/7 — 20 for the violation 
of the Bell inequality is much more stringent than the 
experimenatally realistic rci/7 — 4 that was shown in 
the previous section to be sufficient for producing states 
with (M) > 2. 



VI. CONCLUSIONS 

In conclusion, we have presented a concrete proposal 
for efficient statistical production of multi-qubit GHZ 
states by dispersive measurement in a cQED setup, tak- 
ing into account the realistic conditions of decoherence 
and decay. Our proposal is based on the possibility of 
adjusting the dispersive shifts of individual qubits, which 
effectively modifies the measurement operator and allows 
for the generation of entanglement starting from separa- 
ble input states. Our simulations show that even with 
experimentally achievable values of 2 < rci/7 < 4 it is 
possible to achieve a 1% efficiency in preparing states 
with values of the Bell-Mermin operator exceeding its 
classical bound, (M) > 2. 

By using the global dispersive measurement in the 
same setup, we have also proposed a scheme for imple- 
menting parity measurements on the prepared state. Us- 
ing these measurements, we have studied the sufficient 
conditions for verifying that such states indeed violate 
the Bell-Mermin inequality. We find that a ratio of 
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rci/7 = 20 (essentially identical to the signal-to- noise ra- 
tio) will be sufficient to observe a violation of the Mermin 
bound. While this ratio is larger than currently demon- 
strated, we are optimistic that the present limits on de- 
tector efficiencies in semiconductor amplifiers (1/20 of 
the quantum limit) may be improved by using supercon- 
ducting pre-amplifiers 18]. It would be interesting in the 
future to theoretically explore the possibility to use the 
full detector tomography for violation of the Bell-Mermin 
inequality. This could reveal the necessary conditions 
from the measurement setup for disproving local-hidden 
theories, putting less stringent constraints on experimen- 
tal capabilities. 
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